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ABSTRACT 

We present far-IR ISO observations of three early-type galaxies in the Virgo 
cluster. The data were recorded using the ISOPHOT instrument in both the 
P32 oversampled maps and the P39/37 sparse maps. The maps reach the limit- 
ing sensitivity of the ISOPHOT instrument at 60, 90, and 180/xm. Two of the 
most prominent elliptical galaxies in Virgo - NGC 4472 and NGC 4649 - clearly 
show no emission at far-IR wavelengths at a level of few tens of mJy, but NGC 
4636 is detected at all three wavelengths. We have computed the far-IR emission 
from dust for NGC4472 and NGC4636 under the assumption that dusty outflows 
from evolving red giant stars are continuously supplying dust to the interstellar 
medium and that the grains, once diffused into the interstellar medium, are sput- 
tered away by collisions with ions. While the calculated fluxes are consistent with 
the observed upper limits for NGC4472, the dust emission detected in NGC4636 
support the hypothesis that additional dust has been accreted in a very recent 
(£s few 10 s yrs) merger with a dusty, gas-rich galaxy. 

Subject headings: galaxies: elliptical and lenticular; galaxies: ISM; galaxies, in- 
dividual(NGC 4472, NGC 4636, NGC 4649); infrared: galaxies; infrared: ISM 

1. Introduction 

The interstellar medium (ISM) in elliptical galaxies, although mostly X-ray emitting 
hot gas, contains a significant amount of dust and cold gas of uncertain origin. Most of the 
ISM within the half-light radius is produced by stellar mass loss. Since the mean stellar 
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abundance in typical giant elliptical galaxies is comparable to solar, envelopes ejected from 
evolving giant stars should produce dust just as observed in similar stars in the Milky Way. 
In fact, relatively hot circumstellar dust has been observed in many elliptical galaxies at 
~ 10/ini with both IRAS and ISO satellites (Knapp et al. 1989; Ferrari et al. 2002; Athey 
et al. 2002). This mid-IR emission has a de Vaucouleurs profile similar to the background 
starlight and is therefore consistent with circumstellar dust recently ejected from evolving 
red giant stars. After about 10 4 — 10 5 years, the gas and dust ejected by giant stars orbiting 
in elliptical galaxies expands sufficiently to interact with the ambient hot gas. According to 
currently popular assumptions, as the ejected gas is dynamically disrupted, it is (shocked 
and conductively) heated to the virial temperature (~ 10 7 K) of the hot ISM and the new 
dusty gas merges with its hot environment. When surrounded by gas at ~ 10 7 K, the dust 
is sputtered until it also merges into the hot gas typically on a longer timescale. However, 
during the sputtering time, which varies with galactic radius, the dust receives enough energy 
from electrons and starlight to emit copiously in the mid and far infrared. 
A substantial number of elliptical galaxies have been observed using the ISOPHOT (Lemke 
et al. 1996) instrument on board the ISO (Kessler et al. 1996) satellite. The spectral region 
covered by ISOPHOT, between 60 and 200/iin, together with its mapping capabilities and 
good angular resolution, make the instrument well suited for detecting thermal emission from 
interstellar dust. In particular the extension of its spectral coverage to longer wavelengths, 
relative to the IRAS satellite, provides a better sampling of the far-infrared emission where 
the expected flux density from the cold interstellar dust component peaks. ISOPHOT also 
has a fainter limiting magnitude compared to IRAS and diffraction limited imaging that 
could help resolve and disentangle the various expected sources of emission: central dust 
clouds, circumstellar dust, and truly interstellar dust distributed throughout the galactic 
volume. Here we report on ISOPHOT observations and analyses of three bright nearby 
early-type galaxies that are very similar in their optical parameters but are significantly 
different in their far-IR emission. 



2. Observations and Data Reduction 

ISOPHOT observations of early-type galaxies have been taken in several observing 
modes (Astronomical Observing Template, referred to as AOTs), including the oversam- 
pled maps P32, the P37/38/39 sparse maps and P22 multi-filter photometry data. Here 
we present data from three very bright ellipticals located in the Virgo cluster: NGC4636, 
NGC4649 and NGC4472. 

NGC4636 has been observed using the CI 00 3x3 and the C200 2x2 detector arrays with 
the sparse map on-off mode (AOT P37/39) at 60 and lOO^m (C100) and at 180/wn (C200). 
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Observations were made in each filter with one on-source single point staring mode and one 
off-source exposure. The background position was located ~ 12' north of the target, on a 
blank sky position. Data reduction and calibration were carried out using the PHT Interac- 
tive Analysis tool PIA, version 9.1 (Gabriel et al. 1997) together with the calibration data 
set V4.0 (Laureijs et al. 1998). The reduction included correction for the non-linear response 
of the detectors, readout deglitching, and linear fitting of the signal ramps. After resetting 
of all ramp slopes and subtracting the dark current, the calibration of detector responsivity 
and its changes was performed using measurements of the thermal fine calibration source 
(FCS) on board. The sky subtracted on source signal was then corrected for the fraction of 
the point-spread function included within the field of view of the detectors. Finally, the flux 
densities were extracted applying a color correction assuming that the SED of the galaxy can 
be approximated with a blackbody of temperature 30 K. The statistical errors derived from 
signal processing are about 5-20%, depending on the wavelength range and object brightness, 
but systematic errors due to absolute calibration accuracy are estimated to be 30% (Klaas 
et al. 1998). The derived flux density for NGC4636 in the three filters is shown in Table 1. 
Values are in good agreement with previous IRAS observations. 

NGC4649 and NGC4472 have been observed in the P32 mode in three broadband filters at 
60, 90 and 180/im using the C100 and C200 detectors. The maps were obtained scanning 
the spacecraft Y and Z axes in a grid of 4x4 (NGC4649) and 5x4 (NGC4472) points. The 
focal plane chopper was stepped at intervals of one-third of the detector pixel at each raster 
position, providing a sky sampling in the Y direction of 15" and 30" for the C100 and C200 
detectors. All the maps were centered on the position of the nuclei of the two galaxies cov- 
ering a field of about 9' x 6'. To reduce the P32 data we used a dedicated software package 
(P32TOOLS) developed at MPI Kernphysik in Heidelberg, and supported by the VILSPA 
ISO Data Center (Tuffs & Gabriel 2002). The new routines allow a proper correction for 
transients in PHT32 measurements. The ISOPHOT C200 and especially the C100 detectors 
have a complex non-linear response as a function of illumination history on time scales of 
0.1-100 sec, depending on the absolute flux level as well as the flux changes involved. Since 
no theory exists to predict the detector response behavior, an empirical model has been 
found to apply the transient correction to the data. The severity of the problem depends on 
both the source/background ratio and the dwell time at each pointing direction. Generally, 
ISOPHOT data are reduced in one pass from the raw input data with full time resolution 
to the final calibrated map. To correct for responsivity drift effects, however, it is necessary 
to iterate between a sky map and the input data at full time resolution. This different con- 
cept has been implemented in the dedicated data reduction package. The implementation of 
these new algorithms has been already successfully applied to an extended sample of Virgo 
cluster galaxies by Popescu et al. (2002). Great effort has been devoted to the removal of 
spikes and longer lived tails in the detector responsivity that are often generated by cosmic 
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rays hitting the detector pixels. For the observations presented here, which are faint discrete 
sources, an accurate deglitching procedure becomes critical in determining the instrumental 
sensitivity. Although P32TOOLS applies specific deglitching procedures to remove glitches 
from the data, some residuals are still present in the maps made with the C100 detector. 
Removal of these artifacts was fully achieved for the C200 detector. 

The photometric calibrations of discrete sources performed with the ISO flux scale have been 
tested against integrated flux densities measured by IRAS at 60 and 100/xm bands on Virgo 
cluster galaxies (Tuffs et al. 2002). A very good linear correlation is found between the ISO 
and IRAS fluxes [see Figures 7 a and 7b in Tuffs et al. (2002)]. 

Elliptical galaxies NGC4472 and NGC4649 were not detected at any of the three wavelengths 
observed by ISOPHOT (60, 90, and 180/im). Figure 1 shows the 90/zm maps for these two 
galaxies; the bright source at the edge of the NGC4649 map is its spiral companion NGC4647 
and it is prominent at all three wavelengths. The angular resolution of the oversampled maps 
clearly resolves NGC4647 showing that there is no contribution to the signal associated with 
NGC4649. In the literature, IRAS detections at 60 and 100/im have been erroneously at- 
tributed to NGC4649, simply because the spatial resolution of the IRAS observations, 3'- 
5' at 100/im, is very large and comparable to the optical size of most elliptical galaxies. The 
maps reach the limiting sensitivity of the ISOPHOT instrument at a level of a few tens of 
mJy and the 3a upper limits are presented in Table 1. 



3. A MODEL FOR NGC 4472 AND NGC 4636 

We estimate the FIR emission from dust dispersed in the ISM along the lines of Tsai & 
Mathews (1995, 1996). A full description of the method will be presented elsewhere (Temi 
et al. 2003) In brief, we assume that the grains are ejected from evolving stars and mix with 
the hot ISM in a short time, ~ 10 5 - 10 6 yr (Mathews 1990). Once diffused into the hot gas, 
the grains are sputtered at a rate 

^ = -3.2 x 10~ 14 n p [l + (2 x lOVT) 2 - 5 ]' 1 /mi s" 1 , (1) 
at 

where a is the grain radius (in /im) and n p is the proton density (Draine & Salpeter 1979; 
Tsai & Mathews 1995) 

Let N(r, a)da be the steady state number of grains per cm 3 with radius between a and 
a + da at galactic radius r. Tsai & Mathews (1995) showed that the grain population must 
satisfy 
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where S(r,a)da is the rate at which grains with radii a are expelled from stars per cm 3 . 
Mid-IR (6 - 15/im) ISO observations of a sample of elliptical galaxies, including NGC4472, 
detect strong dust emission that can be understood as emission from circumstellar dust from 
mass- losing AGB stars at or near their globally expected rate (Athey et al. 2002). We assume 
S(r, a) = S (r)a~ s for < a < a max , with s = 3.5 (Mathis et al. 1977), so the coefficient 

Sow = A 36a * p :J 4 ;J 0) 

AnpglO a max 

depends on the rate that the old stellar population ejects mass a* = 4.7 x 10~ 20 s _1 , the 
stellar density p*(r), the initial dust to gas mass ratio scaled to the stellar metal abundance 
5 = (l/150)z(r), and the density of silicate grains, p g = 3.3 gm cm -3 . The solution of 
Equation (2) is then 



N(r, a) 
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a s a< a max . (4) 



The grain temperature is determined by the balance of heating (by both absorption of 
stellar radiation and by electron-grain collisions) and cooling: 

POO 

/ 47rJ*(r, \)Qabs(a, X)-Ka 2 d\ + Ana 2 (l/A)n e (v e E e )T(a) = 4ira 2 a S BTd(r, a) A (Q abs )(T d , a). 
Jo 

(5) 

Since grains experience a relatively high collision frequency with electrons and photons, we 
have not considered stochastic temperature excursions of individual grains (Tsai & Mathews 
1995). Since the 9.7/im silicate feature is seen in emission in a number of elliptical galaxies, 
including NGC4472 (Athey et al. 2002), we assume that all grains have properties similar 
to astronomical silicates (Laor & Draine 1993). Moreover, to simplify the calculations we 
approximate the absorption coefficient with Q a b s pa aip(X) and the Planck-averaged values 
with (Qabs) ~ 1-35 x 10~ 5 Tja [see Draine & Lee (1984)]. The mean intensity of starlight 
is found by integrating over a de Vaucouleurs stellar profile modified with a core of slope 
oc r~ p within r c . At a distance d — 17 Mpc, NGC4472 has an effective radius R e = 8.57 
kpc, total stellar mass M* )t = 7.26 x 10 11 M [assuming M %t /L B = 9.2 (van der Marel 1991)], 
core radius r c = 200 pc and p = 0.9 (Faber et al. 1997). For NGC 4636, assuming again 
d=17 Mpc, R e = 8.32 kpc, M M = 3.80 x 10 11 M , M^ t /L B = 10.74, r c = 260 pc, p = 1.0. 
The collisional heating term in Equation (5) includes a correction r(a) for small grains 
which do not completely stop the electrons (Dwek 1986) and (v e E e ) = y / 32/™^(fcT) 3 / 2 is 
an average over a Maxwellian distribution. The ISM electron density and temperature are 
calculated using the fitting formulae given in Brighenti & Mathews (1997). 

Finally the FIR emissivity is calculated from 

j(r, A) = — / N(r, a)47ra 2 W' 8 Q abs (a, A) x irB(T d (r, a), X)da. (6) 
4tt Jo 



- 6- 



Figure 2 shows the resulting dust temperature profile for NGC 4636, for grains of three 
different sizes. Slightly higher temperatures are predicted for NGC 4472. We find that 
collisional heating dominates over stellar radiation heating for grain sizes a < 0.3 /im and 
a < 1 fim for NGC 4472 and NGC 4636 respectively. 

The computed far-IR spectra calculated for a max = 1 /im and a max = 0.1 /jm are shown 
in Figure 3, together with the ISO detections or upper limits. For NGC 4472 the predicted 
fluxes are consistent with the ISO upper limits provided a max < 1 /im. Thus, no strong 
constraints on the properties of interstellar dust can be derived for this galaxy. On the other 
hand, for NGC4636 the detected fluxes are a factor ~ 50 above the computed values at all 
three observed wavelengths, assuming a max = 1 /im. The shape of the calculated spectral 
energy distribution, however, fits the observations quite well, indicating that the excess 
emission is caused by a larger amount of dust rather than by hotter grains. The calculated 
flux ratios F 90fim / F 60fim = 2.7 and 3.8 (a max = 0.1 and 1 /im); F 180fim / F 90fim = 1.2 and 2.1 
(amax = 0.1 and 1 /im) compare reasonably well with the observed values F 99tira / F m ^ ni ~ 2.6 
and -Fi80//m/^90//m ~ 1-6. The model for NGC 4636 predicts a total mass of interstellar dust 
M dust = 1.7 x 10 5 M© if a max = 1 /im or M dust = 1.7 x 10 4 M for a max = 0.1 /im. We are 
led to speculate that NGC 4636 has accreted dust in a recent merger with a gas-rich galaxy. 
Other explanations appear less palatable. A higher (by a factor > 50) dust production rate 
by a younger stellar population would be inconsistent with the optical colors of NGC4636. 
Moreover, we found that the underluminosity of internally produced dust with respect to 
the observed values cannot be eliminated by varying the grain size distribution parameters 



4. Conclusions 

We have presented ISO observations of three bright elliptical galaxies in the Virgo 
cluster. While NGC 4472 and NGC 4649 emit no far-IR emission at the limiting sensitivity 
of the instrument, NGC 4636 is clearly detected. These observations have been interpreted 
with theoretical models for the far-IR emission from interstellar dust under the assumption 
that the circumstellar dust observed in the mid-IR (Athey et al. 2002) ultimately becomes 
interstellar and is exposed to the mean galactic starlight and to bombardment by thermal 
electrons in the hot gas. Grains receive somewhat more energy from electrons than from 
stellar photons. The calculated far-IR spectrum peaks at A ~ 180 /im, a spectral region not 
accessible to IRAS, but well covered by ISO. 

Although NGC 4472 and NGC 4636 are similar in terms of mass, luminosity, hot gas content, 
and other physical parameters, they show a remarkable discrepancy in their far-IR emission. 
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While the upper limits for NGC 4472 are consistent with emission from dust produced by 
the evolving stellar population, the model for NGC 4636 falls short by ~ 50 in explaining the 
observed fluxes. If the dust mixture would include a significant amount of graphite grains, 
the FIR emission would be at most a factor of few larger (Tsai & Mathews 1996), still much 
too low to explain the observations. This strongly suggests that dust has been accreted by 
NGC4636 in a recent merging event with a gas-rich galaxy. A dust mass larger by a factor of 
about 50-70 is needed in order to match the observed far-IR luminosity, if its relative spatial 
distribution is similar to the one calculated. If the required extra-dust has an external origin, 
however, its distribution may be irregular and patchy. Without spatial information on the 
FIR emission we are not in the position to constrain the spatial location of the external dust. 
Ravindranath et al. (2001) find no compelling evidence for dust extinction in NGC4636 but 
the morphological properties of its stellar core are deviant. An asymmetrical distribution 
of diffuse Ha + [Nil] emission with an irregular velocity field has been detected by Caon, 
Macchetto & Pastoriza (2001) in the central region of NGC 4636. The very disturbed hot 
ISM (Jones et al. 2002) also suggests a recent merging event and this is further supported 
by the apparent failure of the hot gas within ~ R e /2 = 4 kpc to be in hydrostatic equilibrium 
(Brighenti & Mathews 1997). 

The hypothetical merger must be very recent because of the short sputtering time. ^From 
Equation (1) the sputtering time for grains of radius a M (in m/x) is t spu u ~ a^/\da^/ dt\ ~ 
1.2 x lO^n" 1 yrs. Since the gas density in NGC4636 varies as n e ~ O.^Sr^ 36 cm~ 3 
(Brighenti & Mathews 1997), we find t sputt ~ 1-2 x 10 8 a At (rfc pc /i? ei fc pc ) 1 ' 36 yrs. The total mass 
of dust that merged is ~ 10 6 — 10 7 M & and the associated gas would have mass ~ 10 s — 10 9 
M Q . For this gas to have been heated to ~ 10 7 K and for some of the dust to have survived, 
the dusty merger must have occurred less than t sputt) or a few 10 8 yrs ago. The total mass 
of hot gas in NGC4636 within ~ i? e /2 where hydrostatic equilibrium fails, ~2x 10 s M , 
is comparable to that of the added gas, so considerable dynamical disruption is possible. 
Hydrostatic equilibrium would be reestablished in a few dynamical times which is ~ 3 x 10 7 
yrs at r ~ R e /2 = 4 kpc. Dust survival and continued dynamical activity both indicate that 
an important merger may have occurred within the last few 10 8 yrs. It would be interesting 
to seek evidence for the merged stars in the stellar spectrum of NGC4636. 

This study is supported by NASA grant 399-20-01 for which we are very grateful. WGM 
is supported by grants from the NSF and NASA. 
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Fig. I. — Optical image (contours) superimposed to ISOPHOT 90/im map for NGC4649 
(left panel) and NGC4472 (right panel). Both galaxies have not been detected at 60, 90, and 
100/xm. The bright source at the edge of NGC4649 map is its spiral companion NGC4647. 

Fig. 2. — Dust temperature for NGC4636 due to starlight and electron-grain collisions com- 
puted as a function of radial distance. For grains with size a < 0.3 /zm the collisional heating 
is larger than radiative heating. 

Fig. 3. — Computed Far-IR flux from NGC4636 and NGC4472 calculated for a max = 1/xm 
(solid lines) and a max = 0.1//m (dashed lines). The ISO detections are shown as open 
triangles for NGC4472 (upper limits), and open circles for NGC4636. 
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Table 1. ISOPHOT photometry. 



Galaxy 


log (L B ) 


log (L x ) 


A 


Flux (ISO) 


Flux (IRAS) 




(Lb©) 


(erg s" 1 ) 


(urn) 


(mJy) 


(mJy) 


NGC 4649 a 


10.73 


41.21 


60 


<137 c 


900 ± 65 d 








90 


< 85 c 










100 




2310 ± 270 d 








180 


<110 c 




NGC 4472 a 


10.90 


41.66 


60 


<147 c 


± 66 








90 


< 99 c 










100 




± 94 








180 


< 87 c 




NGC 4636 b 


10.55 


41.58 


60 


187 ± 57 


200 ± 34 








90 


491 ± 64 










100 




560 ± 179 








180 


790 ± 71 





a Observed in the P32 oversampled map mode 
b Observed in P37/39 one position sparse map mode 
c No detection, 3 a upper limit flux 

d The aperture photometry includes the spiral companion NGC4647 



